SWAT/02/327 
January 2002 



Mesonic wavefunctions in the 
three-dimensional Gross-Neveu model 



S.J. Hands", J.B. Kogut'' and C.G. Strouthos" 

" Department of Physics, University of Wales Swansea, 
Singleton Park, Swansea, SA2 8PP, U.K. 
^Department of Physics, University of Illinois at Urbana- Champaign, 
Urbana, Illinois 61801-3080, U.S.A. 



Abstract 

We present results from a numerical study of bound state wavefunctions in the (2+1)- 
dimensional Gross-Neveu model with staggered lattice fermions at both zero and nonzero 
temperature. Mesonic channels with varying quantum numbers are identified and analysed. 
In the strongly coupled chirally broken phase at T = the wavefunctions expose effects due 
to varying the interaction strength more effectively than straightforward spectroscopy. In 
the weakly coupled chirally restored phase information on fermion - anti-fermion scattering 
is recovered. In the hot chirally restored phase we find evidence for a screened interaction. 
The T = chirally symmetric phase is most readily distinguished from the symmetric phase 
at high T via the fermion dispersion relation. 
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1 Introduction 

Since the development of lattice gauge theory, many efforts have been made to measure the 
mass spectrum of hadrons. The spatial structure of hadrons, on the other hand, has not 
been as well-studied numerically, although experimentalists have measured quantities like 
the charge radius of the pion [Q]. Numerical investigations of hadronic wavefunctions were 
initiated in 1985 by Velikson and Weingarten and immediately afterwards they were fol- 
lowed by various other groups ^ . Knowing the form factor of the hadrons would allow one 
to calculate their radii and higher moments. Hadronic wavefunctions are also useful in the 
calculation of masses because they can be used as trial wavefunctions in the construction of 
hadronic operators with improved overlap on the ground state at small Euclidean time sepa- 
ration. Furthermore, at high temperature DeTar pointed out the possibility of "confinement" 
at large distance scales ^ . Mesonic wavefunctions at non-zero temperature were investigated 
recently on an anisotropic lattice in the quenched approximation by the QCD-TARO Col- 
laboration 1^. Their results suggest that above there can be low energy excitations in 
the mesonic channels that are metastable bound states i.e. the quark and anti-quark tend 
to stay together at least for Euclidean time scales ~ 1/T. These quasiparticles would be 
characterized by a mass scale given by the location of a peak in the corresponding spectral 
function. This picture is consistent with mean field calculations in the Nambu— Jona-Lasinio 
model 1^ which predict that the mesons associated with the chiral order parameter, namely 
the sigma and the pions are "soft modes" just above T^, i.e. their fluctuations acquire a large 
strength with a small width in the spectral density, which implies the existence of long-lived 
quasiparticles. It has also been shown that the spatial wavefunctions for various mesonic 
channels in the high temperature phase Q are strongly localized, because of the nontrivial 
magnetic interaction of colour currents in the quark-gluon plasma. 

In this paper we present results from a study of mesonic wavefunctions measured in nu- 
merical simulations of the [/(iV/)y-invariant three-dimensional Gross-Neveu model (GNM3) 
with both discrete Z2 and continuous U{1) chiral symmetries at both zero and non-zero tem- 
perature. The model with U{1) chiral symmetry is described by the following semi-bosonized 
Euclidean Lagrangian density: 



We treat ifji, ipi as four-component Dirac spinors and the index i runs over Nj fermion species. 
In the case of a Z2 chiral symmetry the vr fields are set to zero in the Lagrangian. At tree level, 
the fields a and vr have no dynamics; they are truly auxiliary fields. However, they acquire 
dynamical content by dint of quantum effects arising from integrating out the fermions. The 
model is renormalizable in the expansion unlike in the loop expansion |T^. Apart from 
the obvious numerical advantages of working with a relatively simple model in a reduced 
dimensionality there are several other motivations for studying such a model: (i) at T = 
for sufficiently strong coupling > it exhibits spontaneous chiral symmetry breaking 
implying dynamical mass generation for the fermion, the pion field vr being the associated 
Goldstone boson; (ii) the spectrum of excitations contains both baryons and mesons, i.e. 
the elementary fermions / and the composite // states; (iii) the model has an interacting 
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continuum limit at the critical value of the coupling, which has a numerical value g"^ /a ~ 1.0 
in the large- A^j limit if a lattice regularisation is employed 0] ; (iv) numerical simulations of 
the model with chemical potential /i 7^ show qualitatively correct behaviour, unlike QCD 
simulations [|12], 0. 

Let us briefly review the physical content of the model as predicted by the large-A^/ 



approach [110, [Til]. For g"^ > g1 the fermion has a dynamically generated mass Mr equal. 



up to corrections of order 1/iV/, to the scalar field expectation value (a) = g'^ (ipip) . The 
inverse Mj^ defines a correlation length which diverges as {g^ — g^)~'^ with the critical 
index u = 1 + 0{1/Nf). As a result of // loop corrections the a and vr fields acquire non- 
trivial dynamics. In the case of the t/(l)-symmetric model, the pseudoscalar vr couples to 
a Goldstone mode, and its propagator has a massless pole in the chiral limit rriQ 0. For 
both U{1) and Z2 chiral symmetries the scalar a, by contrast, has mass 2Mf in the large- A^/ 
limit, and rather than exhibiting an isolated pole, has a continuum of // states extending all 
the way down to this threshold, implying that if truly bound, its binding energy is 0{1/Nf) 
at best. To our knowledge there have been no analytic calculations of the binding energy 
in this channel. Since all residual interactions are sub leading in 1/Nf, we surmise that all 
other mesons are similarly weakly bound states of massive fermions, and hence effectively 
described by a two-dimensional "non-relativistic quark model". For g"^ < g^ the model is 
chirally symmetric, and hence all states massless as mo — >^ 0. A dimensionful scale is still 
defined, however, by the width of a resonance in // scattering in the scalar channel; this 
diverges as {g1 — g"^)^^ with the same exponent v. In this case we have no equivalent of the 
quark model to assist interpretation of the wavefunctions to be discussed below. 

The large- A^/ approach []T^ also predicts that for g"^ > g'^, chiral symmetry is restored 
as temperature T is raised beyond Tc = M//21n2, or for chemical potential fi > Hc = 
Mf. The phase diagram of GNM3 with various global symmetries at non-zero temperature 
and density has been studied extensively in [12—17]. More specifically, it has been shown 
that the thermally induced phase transition of the Z2-symmetric model belongs to the two- 
dimensional Ising universality class and the f/(l)-symmetric model undergoes a Berezinskii- 
Kosterlitz-Thouless transition [l^, in accordance with the dimensional reduction scenario. 



which predicts that the long-range behaviour at the chiral phase transition is that of the 
(d — l) spin model with the same symmetry, because the IR region of the system is dominated 
by the zero Matsubara mode of the bosonic field. An interesting aspect of GNM3 is that it 
manifests both bulk and thermally-induced chirally-symmetric phases, and one of the goals 
of the present study is to compare and contrast their properties. 

The shape and size of a hadronic state M can be observed through the equal time Bethe- 
Salpeter wavefunction given by: 

vl/(x,t) = I dy(OlV^(y,t)^(y + x,t)lM). (2) 

One can extract from the correlation function C(x, t) which is a convolution of the quark 
propagator Gq and the anti-quark propagator Gq and is given by: 

C(x, t)=J dydy,dy2{0myiMy2)Gq{y, t; yi, 0)rG,-(y + x, t; y^, 0)110) (3) 

The Dirac matrix F selects the appropriate spin and parity quantum numbers for the meson, 
i.e., F = 1 for the spin zero scalar (S) and F = 75 for the spin zero pseudoscalar (PS). $(x) 
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is an input trial wavefunction which is used as a source for the construction of quark and 
anti-quark propagators. At large t the contribution from the ground state dominates and 

C(x, t) ~ exp(-mHt)^(x), (4) 

where tuh is the hadron mass. As already mentioned in |^ these wavefunctions are minimal 
Fock space wavefunctions, because they do not overlap onto states for which the quark anti- 
quark world line has kinks crossing time t. Therefore, their use to calculate phenomenological 
numbers is an uncontrolled approximation. However, as we shall see, they are convenient 
tools in lattice simulations to study the binding in the various mesonic channels. Due to its 
simplicity, the GNM3 can be studied at both T = and T 7^ on lattices with relatively 
large Lt, which is the main difficulty for QCD simulations. In the next section we will discuss 
numerical results of mesonic masses and wavefunctions on large volumes in both the scalar 
(S) and pseudoscalar (PS) channels. In our study we only measured the connected parts 
of the correlators. The noisy disconnected diagrams were neglected and we will discuss the 
implication of this in the next section. We will also discuss the dependence of our results 
on the quark propagator source, on the lattice spatial extent L^, the chiral symmetry group 
{Z2 vs. U{1)) of GNM3, the coupling f3 = l/g"^ and the number of fermion species in the 
model Nf. We will then present results in the T = and T 7^ symmetric phases. The 
sigma and the pion are represented in the semi-bosonised GNM3 action by bosonic auxiliary 
fields and hence the correlation functions in these channels including disconnected diagrams 
can be measured with relatively high statistics. Unfortunately, the same technique cannot 
be applied to the wavefunctions which are point-split quark four-point functions. In another 
project, which in a sense is complementary to this one we are studying mesonic spectral 
functions including those of the auxiliary fields using the Maximum Entropy Method [|1^. 



2 Simulations 

The fermionic part of the lattice action we have used for the semi-bosonized GNM3 with 
U{1) chiral symmetry is given by JT^ 

Sfer = Xi{^)MijxyXj{y) 

^ / 1 \ 

= II IIXi(a;)-^^i/Xi(2/) + ^IIXi(a;)Xi(a;)[I](x(x) + ze(a;)^7r(£)] , (5) 

i=l \^,y ° X (x,x) {x,x) I 

where Xi and Xi are Grassmann-valued staggered fermion fields defined on the lattice sites, 
the auxiliary fields a and vr are defined on the dual lattice sites, and the symbol (x, x) denotes 
the set of 8 dual lattice sites x surrounding the direct lattice site x. The fermion kinetic 
operator M. is given by 

= ^ H ^'^(^) ~ ^y,x-A + "^o^xj/, (6) 

where ri^{x) are the Kawamoto-Smit phases (— 1)^°"' '"^"-1^ and the symbol e{x) denotes 
the alternating phase (^— i)^o+^i+^'2. The auxiliary fields a and tt are weighted in the path 
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integral by an additional factor corresponding to 



Sau. = ^,j:W'ii) + 7r'{x)]. (7) 

The simulations were performed by using the standard hybrid Monte Carlo algorithm in 
which complex bosonic pseudofermion fields $ are updated using the action $^(M^M)~^$. 
Unless stated otherwise, the bare fermion mass mo was set to zero. According to the dis- 
cussion in |12], simulation of staggered fermions describes Nf = AN continuum species; 
the full symmetry of the lattice model in the continuum limit, however, is U{Nf/2)y ® 
U{Nf/2)v ® U{1) rather than U{Nf)y ® U{1). At non-zero lattice spacing the symmetry 
group is smaller: U {Nf/4)v ® U{Nf/4)v ® ^(1)- In the Z2-symmetric model the vr fields are 



switched off and M becomes real. In this case according to the discusion in |]TT| simulation 
of A^ staggered fermions describes Nf = 2N continuum species. Details about the algorithm 
and how we optimized its performance can be found in [jrT], |T^. 

Using point sources we calculated the zero momentum fermion correlator at different 
values of the coupling /3 and fitted to the function 

Cf{t) = Af[e-^'f' - (-l)*e-*'^^(^*-*)]. (8) 

The mesonic correlators are given by: 

CM{t)= J2 <f(xi)$(x2)W^M(x)G(x,t;xi,0)Gt(x,t;x2,0), (9) 

X,Xl,X2 

where Wm{'^) is a staggered fermion phase factor which picks out a channel with particular 
symmetry properties i.e. PVm(x) = e(x) for the S channel and PVm(x) = 1 for the PS 
channel. The function $(x) is either a point source (5x,(o,o) or a staggered fermion wall source 



J2m,n=o ^x,(2m,2n) EH- Ih all the simulatious we used point sinks. These correlators were 



fitted to a function C^/ (t) given by 

CM{t) = Ale-""' + e-'"^^'-*)] + i(-l)*[e-™* + e-"'^^'~% (10) 

The first square bracket represents the "direct" signal with mass m and the second an 
"alternating" signal with mass fh. Just as in four dimensions, composite operators made 
from staggered fermion fields project onto more than one set of continuum quantum numbers. 
To gain more insight one must transform to a basis with explicit spinor and flavor indices 
|20| ; the resulting bilinears g(r T)q in both S and PS channels, where T acts on a two- 



component flavor space, are summarised in Table |I|. Also shown are the spin/parity quantum 



Table 1: Spin/flavor assignments of the fermion bilinears studied in the paper 
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numbers . In the direct channels these are the same as their four- dimensional equivalents. 
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but the alternating channels are very different. In the PS channel the alternating state is 
an anti-symmetric tensor, corresponding to spin-1, whereas in the S channel the bilinear 
transforms trivially under rotations in the 12 plane and is hence spin-0. The parity states 
are assigned by considering transformation properties under 2+1 dimensional lattice parity: 

X = {xo,Xi,X2) 

a{x) I— s> (t(x') 

In continuum notation this corresponds to determining whether F commutes or anti-commutes 
with 7i75. It is important to note that unlike in 3+1 dimensions, the alternating channels 
do not contain "parity partners" of the direct channels, or indeed necessarily even decribe 
states of the same spin. 

In a similar way the lattice wavefunctions are given by: 

^(x)= J2 $(xi)<l'(x2)W^M(x)G'(x + x',t;xi,0)G't(x',t;x2,0). (12) 

x',Xl,X2 

In the staggered formulation the easiest way to construct is by computing it for spatial sep- 
arations that are multiples of 2a. To facilitate comparison between the various wavefunctions 
we normalize them to unity at zero separation. 



X = (Xo, 1 - Xi,X2) 
7i(x) ^■ 



(11) 



2.1 Zero Temperature Broken Phase 

In this section we discuss the numerical results for the masses and the wavefunctions ex- 
tracted from the S and PS channels of GNM3 at T = 0. We generated statistics that vary 
from 700 to 2000 configurations, depending on how close to the critical point is the coupling 
for each simulation. The average molecular dynamics trajectory length f ~ 1.2. 

First, we measured the meson correlators Cuit) in the S and PS channels and the fermion 
correlator Cf{t) at different values of the coupling (3 from simulations of the Z2-symmetric 
model. An interesting observation is that the direct amplitude A in the S correlator is tiny 
compared to the alternating signal A, and appears to decrease rapidly as we decrease the 
coupling /?, whereas the PS correlator is dominated by the direct signal. For example, in the 
case of the S correlator for [3 = 0.75 the fit parameters are A = —0.00010(5), A = 250(10), 
m = 0.34(1), Mf = 0.175(1) and for p = 0.80 they are A = -0.05(1), A = 235(5), 
m = 0.172(1), Mf = 0.0893(5). The bulk critical coupling is Pc = 0.85(1). We infer 
that the alternating signal corresponds to a weakly bound state since its mass is slightly less 
than 2Mf. It appears that the flavor- singlet scalar channel Sdir is completely dominated by 
disconnected diagrams in this phase, and can only be studied using the auxihary a field [pT| . 
Because of the smallness of A in the S correlator {A in the PS correlator) for /3 ^ /3c — 0.85, 
and the statistical fluctuations in our data we do not quote values for their masses; they 
did, however, appear very light, perhaps even consistent with zero, perhaps suggesting their 
origin is due to the proximity of the chirally symmetric phase at /3 ^ 0.85. 

Even in the case of continuous chiral symmetry, the meson extracted from the connected 



contribution in the PSdir channel has a mass almost twice the fermion mass O, p^Sl because 
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it is not the physical pion: the Goldstone mechanism in GNM3 is fundamentally different 
from that in QCD. In this model the disconnected diagrams are responsible for making the 
pion light; numerically this can be accessed by studying correlators of the auxiliary vr field 
0,111. 

In Fig. |1| we present results for the PS channel wavefunctions of the Z2 model with 
Nf = 4 on a 32^ x 96 lattice and coupling (3 = 0.55. The wavefunctions were extracted 
from different timeslices {t = 11,21,31,41) and were constructed from quark propagators 
with both wall and point sources. It is evident that the wavefunctions extracted from wall 
sources stabilize at relatively small Euclidean times. It is expected that the wall source 
projects onto the ground state much more effectively than the point source, because the 
wavefunction is very broad and it does not go to zero at the edge of the lattice, i.e. it 
resembles a wall wavefunction. The other advantage of the wall sources wavefunction is that 
it is constructed by convolving quarks with p = 0, whereas the wavefunction extracted from 
point sources has more contamination from excited states at small times. As mentioned in 
the previous paragraph this weakly bound state has a mass which is very close to twice the 
fermion mass; in this case Mps = 1.085(3) and Mf = 0.550(1). Because of the superior 
saturation properties, the various correlators that we will discuss in the next few paragraphs 
were constructed from wall sources. 

In Fig. ||we present PSdir and Salt wavefunctions at coupling (3 = 0.55 and lattice sizes 
32^ X 96 and 48^ x 96. The wavefunctions fill the lattice even when Lg = 48 and it is obvious 
that finite size effects are still large at the boundary. In Fig. | we plot the Salt and PSait 
wavefunctions in the broken phase for (3 = 0.45, 0.55, 0.65 with the horizontal axis rescaled by 
the fermion mass Mf which is an inverse correlation length. For small distances (in physical 
units xMf ^ 3) the three PSdir wavefunctions coincide, implying that in this regime the 
wavefunction shape is that of the continuum limit theory. We interpret the deviations at 
large distances as strong finite size effects; as implied by Fig. |^, the meson fills the lattice 
at these couplings. The Salt wavefunctions which are broader than the PSait wavefunctions 
deviate slightly from each other even near the origin, indicative of stronger curvature in this 
region. 

One clear effect is that the // pair is much more closely bound in the PS channel than 
in the S. We can explain this by adapting the well-known argument from the non-relativistic 
quark model that a fermion - anti-fermion system bound in an s-wave is necessarily pseu- 
doscalar due to the product of intrinsic parities of the // pair being negative. Some care 
must be taken in 2+1 dimensions; whereas the L = spatial wavefunction is axially sym- 
metric and hence positive parity, all higher angular momentum states, which have the form 
f{\r\)e^'^^, where /(r) oc r for small r and for non-singular potentials, contain both pos- 
itive and negative parity components (unlike in 3+1 dimensions where parity is given by 
(— 1)"^), implying that higher spin states fall into degenerate parity doublets if parity is a 
good quantum number Q. The conclusion remains, however, that the only allowed s-wave 

"'^ Intrinsic parity must also be defined with care: recall that in 3+1 dimensions it is related to the spinor 
eigenstates of 70. Since the continuum formulation of staggered fermions in 2+1 dimensions naturally results 
in 4-component spinors, we can effectively regard the parity transformation (|l^) as a spatial inversion in 3 
dimensions followed by a rotation by tt radians in the 23 plane; since rotations are continuous the / and / 
spinor eigenstates maintain their respective intrinsic parities under this sequence of transformations. 
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state is the 0~, which naturally accounts for why the FS^ir wavefunction has the smallest 
extent. The behaviour of states with L > as a function of |r| might possibly explain the 
noticeable kink in the S wavefunction at the smallest non-zero spatial separation; ideally 
data nearer the continuum limit would be needed to explore this more thoroughly. 

In Fig. ^ we present PSdir and Salt wavefunctions measured in simulations with different 
numbers of fermion species in the system, i.e., Nf = 4, 12. The coupling is set to (3 = 0.55 and 
the wavefunctions shown in this figure were extracted from the t = 40 timeslice. In order to 
compare the two cases and probe 1/Nf effects we rescaled the horizontal axis by multiplying 
the distance with the fermion mass Mf. It is clear from the figure that as we switch on the 
interaction (i.e. decrease Nf) the binding increases and the wavefunction shrinks. The ratio 
Mps/{2Mf) is 0.986(5) for Nf = A and 0.999(1) for Nf = 12. These results imply that the 
wavefunctions in this model are far more sensitive to 1/iV/ effects than the masses or even 
the critical indices. In this model the bulk critical indices for finite Nf are very close to the 
large- A^/ prediction and 1/A^/ effects have been very difficult to measure even for relatively 
small Nf [ID. 

In Fig. ^ we present wavefunctions measured in the PS^j^ and Salt channels of GNM3 
with either a Z2 or a U{1) chiral symmetry. In both cases the number of fermion species is 
Nf = 4 and the coupling is (3 = 0.55. In the U{1) model we fixed the vacuum by introducing 
a small fermion bare mass mo = 0.01 so that (a) > 0. The masses in the U{1) model are: 
1.133(3) for the S channel, 1.123(3) for the PS channel and 0.567(2) for the fermion. In the 
Z2 model the masses are 1.085(5) for the S channel 1.083(5) for the PS channel and 0.550(2) 
for the fermion. These mesons have very broad but quite distinct wavefunctions, whereas 
their masses are almost the same and are very close to twice the fermion mass. This is 
further evidence that the wavefunctions in this model are more sensitive to binding effects 
than bound state masses. The difference between PS and S is larger in the U{1) model 
than in the Z2 model: the PS is more tightly bound, due to the additional attractive force 
resulting from vr exchange, which since the vr is lighter moreover also results in a longer-range 
// potential. The kink in the S wavefunction is also more pronounced in this case. Although 
the various meson masses measured in our simulations are in accordance with the large- A^j 
prediction slightly less than 2Mj, we believe that an accurate measurement of the binding 
energy can only be extracted from simulations on lattices which are bigger than the meson 
sizes, otherwise finite size effects can severely distort the result. 



2.2 Zero Temperature Symmetric Phase 

We now consider the wavefunctions in the chirally symmetric phase. As discussed in the 
Introduction, here there is no dynamically generated mass scale, and no readily identifiable 
correlation length, although evidence for a resonance width fi scaling approximately linearly 
in ((yf^ — g^) was reported in In Fig. || we show the PS^j^ wavefunctions of the Z2 model 
in both broken {(3 < (3c) and symmetric {(3 > (3c) phases. The main observation is that the 
wavefunction becomes broader, i.e. the size of the meson increases, as (3 increases, but that 
some // attraction persists beyond the critical point (3c = 0.85(1). Counterintuitively, the 
wavefunction size in lattice units continues to increase for (3 > (3c, whereas the correlation 
length should decrease as we recede from the critical point. This may indicate that there is no 



8 



true bound state in the channel, but merely a positive spatial correlation due to interparticle 
attraction, which grows weaker as /5 — > oo. Before taking this interpretation too seriously, 
however, it should be remarked that taken on their own the data of Fig. |^ would not enable 
the critical coupling to be identified. 

In Fig. we plot data extracted from a simulation on a 32'^ lattice at coupling (3 = 1.0 
in the symmetric phase in both S and PS channels. We generated 10, 000 configurations. 
A simple analysis of the correlators Cm(^) implies that both are massless. Compared with 
the broken phase, however, the saturation is much improved, permitting a more complete 
analysis, which ultimately permits extraction of both direct and alternating signals from 
each channel. It is possible that in this case the wavefunction is yielding information on 
particle scattering rather than on an isolated bound state - this may explain the improved 
saturation, since there is no longer a requirement for the signal from excited states to die 
away as t increases before the shape of "^{x) stabilises. 

We disentangled "^dirix) and "^aitix) by measuring \l'(a;,t), which is given by 

^(X, t) = A^dirix) + i-iyA^altix), (13) 

at large even and odd t. The results presented in Fig. |^ were extracted from the values of 
at t = 14, 15. A striking feature of the plot is that |\E'p5| + |\1'5| appears independent of x. In 
fact, the meson amplitudes are of opposite signs, implying that the sum of factors Wps + Ws 
(see (p!2D) is of the form (1 — e(x)) and hence projects onto odd sink lattice sites only. Now, 
in a chirally symmetric vacuum, meson states built from non-interacting fermions (implying 
a spatially uniform ^f) have contributions of precisely this form GoeGl^; quantum corrections 
to |\E'ps'| + l^'sl can therefore only arise from exchange of an even number of auxiliary bosons, 
and hence are suppressed by 0(1/A^|) at least. As already observed the PS wavefunctions 
expand as we increase f3, with PS^j^ niore tightly bound that PSa/t. The similarity in shape 
may indicate that according to the assignments of Tab. ^ both are permitted to scatter 
in the s-wave. In contrast to the PS case, there is a qualitative difference in S channel 
wavefunctions between broken and symmetric phases. The wavefunction in both S^j^ and 
Salt increases away from the origin, implying that close approach of the // is disfavoured. 
This is consistent with the requirement of positive parity states to have non-zero orbital 
angular momentum. 

2.3 Nonzero Temperature 

Next we discuss the effects of non-zero temperature T. We have studied lattices with tem- 
poral extent = 16, for which the critical point for the Z2 symmetric model is estimated 
to be /3j = 0.790(5) [^, and the high temperature chirally restored phase can thus be 
studied for /3j < /5 < P^'^''''. A particular virtue of GNM3 is that, unlike unquenched QCD, it 
permits a study of hot dynamics with current resources with Lt sufficiently large to permit 
the measurement of thermal masses from correlators in Euclidean time. We have extracted 
the thermal fermion mass Mf by fits to 

Cf{t) = A{1 - i-lY) sinh(Mf (t - Lt/2)) + B{1 + (-1)*) cosh(Mf (t - Lt/2)). (14) 
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We performed simulations on 32^ x 16 and 96^ x 16 lattices and we generated approximately 
20, 000 configurations on the smaller lattice and 35, 000 on the bigger one. However, our 
data are consistent with the vanishing of Cf{t) for even t, corresponding to 5 = in (]T^), 
signalling a manifest chiral symmetry. The results for Mf^ are given in Tab. |^. We observe 



Table 2: Thermal fermion masses in the hot phase 



p 


32^ X 16 


96^ X 16 


0.79 




0.048(1) 


0.80 


0.044(1) 


0.041(1) 


0.82 




0.034(1) 


0.85 


0.026(1) 





that Mf is reasonably insensitive to the spatial volume of the lattice. The masses, though 
small, are considerably larger than those extracted from the bulk symmetric phase on a 32^ 
lattice, where typical values found were Mf = 0.006(1) at /? = 1.0 and 0.0015(3) at /3 = 2.0, 
ie. virtually massless. The fermion mass thus provides a criterion from distinguishing the 
two chirally symmetric phases. In thermal field theory the typical scale for thermal masses 
is given by m'^^ ~ g(T)T where g(T) is the coupling strength associated with a vertex in a 



perturbative calculation. For GNM3 the coupling associated with the UV fixed point is ||TT 



(15) 



For asymptotically high T we thus have g{T) oc T and hence limT-^oo ^/ const. For 
finite T we still expect Mf'/T to increase more slowly than T, implying that for fixed 
Mj'^a should decrease as (3 —>■ (3^^^, as suggested by the data of Tab. 0. In this limit the 
large- A^j expansion suggests that GNM3 becomes an ideal fermi gas [P^. 

We have also examined fermion correlation functions at non-vanishing momenta k = 
27rn/Ls, n = 0,1,2,... and used the energies E{k) extracted from (113) to map out the 



fermion dispersion relation. Results from both hot and bulk chirally symmetric phases are 
shown in Fig. 0, togther with fits to the lattice form 



E{k) = A smh-'{^smH + Mj). (16) 

The fits yield masses consistent with those extracted by fitting directly at zero momentum, 
and in both phases yield values of A, related to the renormalised speed of light, very close 
to 1, implying that the principal physical effect of the hot medium is to generate a non- 
zero thermal mass. A similar observation has been made in quenched QCD 0. The result 
y4 ~ 1 supports our identification of a non-zero screening mass, or "pseudogap", in studies 
of the [/(l)-symmetric model at T > [|r^, where large volume effects precluded using L 



t 

sufficiently large to extract E{k). Modifications of the low momentum part of the dispersion 
relation which is expected to have two branches at T > |^ corresponding to two kinds of 



quasiparticle excitations, the fermion and the hole (known as the plasmino), are not visible 
in our simulations. Much larger lattices are needed to study the structure of the dispersion 
relation at the required low momenta. 
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As in the bulk symmetric phase, it is necessary to have accurate measurements of meson 
masses before wavefunctions can be extracted. On 16 x 32^ at /? = 0.80 we found a mass 
0.075(5) in the PSdir channel, which decreased to 0.038(3) at j3 = 0.85, and masses close to 
if not consistent with zero for PSaZj. In the S channel stable signals were found in both Salt 
(0.090(2)) and Sdir (0.047(2)) at /3 = 0.80, decreasing to 0.051(5) and 0.037(5) hy f3 = 0.85. 
As in the bulk symmetric phase, saturation was found to be quite satisfactory, which may 
again be a sign that there is no bound state pole. The resulting wavefunctions in PS^j^ and 
PSait channels are shown in Fig. |[ The general shapes are very similar to those seen in both 
chirally broken and bulk symmetric phases, with the width of the wavefunction increasing 
monotonically with f3. In principle the x axis should be rescaled with the correlation length 
given by MJ^{T = 0); however, we have not done so since P = 0.85 is extremely close 
to the bulk critical point where finite volume corrections are large, making an estimate of 
Mf{T = 0) by any other means than naively extrapolating from smaller f3 using i/ ~ 1 
impracticable. Instead, we note that the physical scale of the structure in Fig. ^ must shrink 
to zero as /5 — > implying that the region where // interactions are important becomes 

pointlike as T ^ oo. This could perhaps be explained in terms of Debye screening of the 
// interaction due to thermal effects. Once again, we are unable to determine the issue of 
whether the mesons remain bound states for T > T^. without a detailed study of the spectral 
function - the behaviour of \E'(x) for large x is subject to large finite volume effects in this 
coupling regime. 

3 Summary and Outlook 

We have applied a technique originally developed for studying mesons in QCD, namely 
analysis of spatial wavefunctions, in a simpler model which manifests a variety of interesting 
behaviours as both coupling strength and/or temperature are varied. We have successfully 
identified mesonic bound states in the chirally broken phase, and exposed the effects of 
changing the strength of the interaction, either by altering the number of species Nf or the 
chiral symmetry group, much more cleanly than permitted by a study of the spectrum alone. 
Our results are reminiscent of the well-known property of the variational method in quantum 
mechanics that energy eigenvalues are better determined, i.e. less sensitive to small changes 
in the variational parameter, than the wavefunctions themselves. The physical picture is 
a two-dimensional analogue of the quark model, with the s-wave nature of the PS state 
clearly distinguishable. Since the states are weakly bound, they have a larger spatial extent 
than the corresponding QCD mesons, and hence are more prone to finite volume corrections. 
We have also studied the wavefunctions in two distinct chirally symmetric phases, a T = 
phase for (3 > and a hot phase for /5j < /3 < Because of the difficulties in 

reliably extracting a correlation length in either case, our study here has necessarily been 
more exploratory. Nonetheless we have seen evidence for a clear channel dependence in 
the spatial correlation of propagating // pairs, which may now be due to a continuum of 
scattering states rather than isolated bound states. Interestingly, the clearest distinction 
between the two symmetric phases comes from a study of the fermion dispersion relation, 
which yields a non-vanishing thermal mass for the T > study. If we use our expectations 
from the large- A^j approach, then the wavefunctions also provide evidence for screening in 
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the hot phase. 

The difficulties we have faced with finite volume effects are intimately related to the 
details of the spectrum of the model. In a related study we are analysing the spectral 
functions of the model using Maximum Entropy techniques. One advantage is that this 
approach permits analysis of the auxiliary fields, corresponding in lattice QCD parlance to 
the inclusion of disconnected quark line diagrams as well. It is hoped that this complementary 
approach will also help to resolve the issue of whether there are bound states in the chirally 
symmetric phases. 

Finally, another feature of GNM3 is that it permits simulation with a non-zero chemical 
potential /i, making studies of degenerate "quark matter" possible [1^. We are currently 
extending our wavefunction studies to ^ ^ 0, where a sharp Fermi surface manifests itself 
by a \l/(a;) which fiuctuates in sign, demonstrating so-called Friedel oscillations p4 
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Figure 1: PS wavefunctions at /? = 0.55. Filled symbols correspond to wall source; t increases 
from top to bottom. Empty symbols correspond to point source; t increases from bottom to 
top. 
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Figure 2: S and PS wavefunctions at (3 — 0.55 and Lg — 32, 48. 



14 



1.1 



1.0 txA 



0.9 



0.8 



0.7 



0.6 



• • PS P=0.45 

O OS M-45 

♦ ♦ PS p=0.55 
O OS M-55 
A ▲ PS |3=0-65 
A ASP=0.65 



* O 



O. ^ O 



* * « 



■D- 



• » ■ 



10 



15 



Figure 3: S and PS wavefunctions at (5 — 0.45, 0.55, 0.65; the horizontal axis is rescaled by 
the fermion mass Mf. 
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Figure 4: S and PS wavefunctions for Nf — A and A^/ = 12; the horizontal axis is rescaled 
by Mf. 
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Figure 5: S and PS wavefunctions from the C/(l) and Z2 GNM3. 
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Figure 6: PS^j,. wavefunctions in the broken phase {j3 = 0.45, 55, 65, 0.75) and in the sym- 
metric phase (/3 = 1.0, 2.0); /3 increases from the bottom to the top of the figure. 
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Figure 7: S and PS wavefunctions in the zero temperature symmetric phase. 
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Figure 8: PS wavefunctions in the high T symmetric phase. 
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